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Introduction 

C0 2 is the only gas which produces appreciable infrared 
cooling of the Martian troposphere (Goody and Belton, 
1967 ; Crisp, 1990 ; Lindner, 1991 ; Savijarvi, 1991 ; Hour- 
din, 1992); hence, the transmission on Mars of the C0 2 
15 pm band has been well studied (Gal’tsev and Osipov, 
1979; Crisp, 1990; Lindner et at., 1990a,b; Savijarvi, 
1991 ; Hourdin, 1992). Airborne dust also produces 
appreciable cooling, even in periods other than global dust 
storms (e.g. Lindner, 1985, 1993). However, virtually all 
previous studies of Martian radiative cooling have ignored 
the presence of dust when calculating C0 2 cooling rates. 
The effect of dust-C0 2 spectral overlap (hereafter abbrevi- 
ated as dust-C0 2 interaction) was studied by Kondrat’ev 
and Moscalenko (1975) and Egan et al. (1980) for the 
effect on line profiles, by Pallman (1983) for the near- 
surface boundary layer, by Pollack et al. (1979) for invert- 
ing data, and by Habcrle et al. (1982) for circulation 
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models. However, none of these studies published or dis- 
cussed dust-C0 2 interaction. The importance of both C0 2 
and dust in Martian radiative transfer requires a model 
capable of accurately handling the scattering by dust and 
the complex wavelength structure of C0 2 . To simplify 
the problem, most prior work has studied C0 2 and dust 
radiative transfer separately and avoided their over- 
lapping opacities. Here I examine the errors that can result 
from such a treatment. 


Model 

The discrete ordinate method of Stamnes et al. (1988) is 
used to treat the scattering, absorption and emission of 
monochromatic radiation through the Martian atmo- 
sphere. The exponential sum program of Evans et al. 
(1980) converts the banded wavelength structure of C0 2 to 
allow for monochromatic treatment (Freeman and Liou, 
1979 ; Lindner, 1985 ; Lindner et al., 1990a,b). The 15 pm 
band is broken into 16 subintervals to properly account 
for the multiple scattering (Lindner et al., 1990a,b ; Laris 
and Oinas, 1991). Transmission functions for the 15.0 pm 
band of C0 2 are taken from the Iine-by-Iinc model results 
of Gal’tsev and Osipov (1 979). The transmission function, 
Tr, as a function of temperature T, pressure P, and C0 2 
column abundance U is extrapolated from the Gal’tsev 
and Osipov results (subscript G) to temperatures below 
200K by: 

TriT, P,U) = 1 — {(1 — 7>7200K, P, U)) (77200K) 0 ‘}. 

(I) 

The exponential 0.8 is found when the temperature depen- 
dencies for the Pollack et al. (1981) transmission functions 
are recast in this form. Using a modified version of the 
FASCOD transmission model (Clough et at., 1986), the 
accuracy of the transmission functions in equation (1) 
arc confirmed for the range in temperature and pressure 
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present in the atmosphere to 40 km altitude (Lindner et 
al., 1990a,b; Pollack et al., 1990). Local thermodynamic 
equilibrium is assumed (Gierasch and Goody, 1967 ; 
Hourdin, 1992). 

Dust opacities vary from 0.2 to 1.0 for conditions other 
than global dust storms (Pollack et al., 1979). However, 
dust opacities over winter polar latitudes may be slightly 
less (e.g. Lindner, 1990). A Gaussian profile describes the 
vertical distribution of dust, being well mixed to 20 km 
altitude for conditions other than global dust storms 
(Anderson and Leovy, 1978; Zurek, 1982; Korablev et 
al., 1993). The wavelength dependence of the dust opacity 
and the dust single-scattering albedo are given by Toon et 
al. (1977). Scattering of radiation by dust is represented 
by the Haze-L phase function (Toon et al., 1977). Com- 
putational difficulties which accompany highly asym- 
metric phase functions are removed with the Delta-L 
method (Wiscombe, 1977). The emissivity of airborne dust 
is high and has been calculated as a function of wavelength 
from theory and observations (Toon et al., 1977 ; Simpson 
et al., 1981). C 

— ^ Atmospheric composition is taken as 95% gQ 2 (Owen 

et al., 1977). Atmospheric composition may have been 
quite different in past epochs, with C0 2 being perhaps a 
minor constituent (e.g. Lindner and Jakosky, 1985), but 
this study focuses on the present epoch. Season-dependent 
Ct\ abundances (Hess et al., 1980) are corrected for 
cir«. ilation-induced pressure gradient (Haberle et al., 
1979) and elevation (Lindal et al., 1979). The surface 
pressure is 8 mbar at 57°N latitude in late winter. Atmo- 
spheric properties are zonally averaged and are assumed 
azimuthally independent. The region from the surface to 
40 km altitude is broken into 20 layers, each 2 km thick, 
to account for vertical inhomogeneity. The temperature 
profile rises linearly with altitude from 150K at the surface 
to 160K at 10 km, and then falls linearly with altitude to 
130K and 40 km, typical for winter polar conditions (Lin- 
dal et al., 1979 ; Kieffer, 1979 ; Martin, 1984). The albedo 
of the polar cap is assumed to be zero (Kieffer, 1970; 
Smythe, 1975; Wiscombe and Warren, 1980). An average 
ice emissivity of 0.9 is adopted, with the wavelength depen- 
dence given in Ditteon and Kieffer (1979) and Hunt et al. 
(1980). 


Results and discussion 

Total atmospheric cooling rates owing to thermal emis- 
sion by a vertical optical depth of dust, x r , of 0.2, and the 
contribution by various wavelength intervals are pre- 
sented in Fig. 1 for winter polar conditions. The maximum 
cooling owing to dust below 30 km altitude occurs in the 
12- 1 8 pm wavelength interval. The strong dust cooling in 
the 12-18 pm range is owing to the cold temperature, to 
the high optical depth of dust at 12-18 /mi (Toon et al., 
1977), and to the secondary maximum in dust emissivity 
at 12-18 pm (Toon et al., 1977). Hence, COj impedes the 
ability of dust to cool the atmosphere because C0 2 absorbs 
strongly in the 12-18 pm region where dust cools most 
effectively. 

The presence of dust will affect C0 2 radiative transfer 



Dust Cooling Rate (K/Day) 


i Fig. 1. Total infrared cooling rate by dust at winter polar lati- 
: tudes (57°N latitude, L, = 343°) as a function of altitude for 
j dust vertical optical depth, of 0.2, and the contribution from 


each wavelength interval from 5 to 40 pm (wavelength range of 
each bin is labeled in pm) 


by increasing the optical depth, and absorbing light that 
is emitted by C0 2 . Dust is particularly important in the 
radiative cooling in the line wings of C0 2 . Line wings are 
important to cooling in that the strongly absorbing line 
center is often optically thick, which renders the line center 
ineffective at cooling (e.g. Lindner, 1993). Optically thin 
line wings cool effectively, in that photons emitted in the 
line wings can escape and cool the atmosphere. The pres- 
ence of dust affects the line wings by absorbing emitted 
photons that would otherwise escape. 

Hence, dust and C0 2 hinder each other’s attempts to 
cool. Thus, simply adding dust cooling to C0 2 cooling is 
not accurate in those wavelength bands where both are 
important. This competition affects C0 2 cooling more 
than dust cooling. Dust cooling occurs over the entire 
infrared spectrum, which is mostly free of C0 2 absorption. 
However, C0 2 absorption and emission must always com- 
pete with the absorption by dust. Figure 2 shows that the 
cooling is overestimated up to 30% if dust and C0 2 
cooling rates are computed separately and then summed. 
The effect is strongest near the surface where the optical 
depth is largest. C0 2 cooling is optically thick, making it 
very difficult for emitted photons to escape in the lowest 
10 km. Some of the photons emitted by C0 2 line wings 
that would normally escape and cool are being absorbed 
by dust, which results in a lower cooling rate. C0 2 also 
absorbs some of the photons emitted by dust that would 
normally escape and cool. 

The error induced by ignoring dust-C0 2 interaction is 
increased at higher dust loading (Fig. 3). Dust is now even 
more capable of choking off the escape of emitted photons, 
not only those of C0 2 , but also photons emitted by dust 
itself. Dust-C0 2 interaction is highest near the surface for 
all three cases. The relative importance of the line wings 
for cooling is greater at lower altitudes. Since dust has its 
strongest effect on the line wings, dust-C0 2 interaction is 
strongest near the surface. Note that the errors do not 
change as much when dust opacity is changed from 
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15 Micron Region Cooling Rate (K/Day) 

Fig. 2. Total 15 /on region (12-18 pm) cooling rates as a function 
of altitude for a pure COj atmosphere (labeled “gas only”) and 
for a pure dust (t, = 0.2) atmosphere (labeled “dust only") at 
winter polar latitudes. Also shown are the cooling rates for a 
simple sum of the C0 2 only and dust only cooling rates (labeled 
‘.‘gas only + dust only"), and the cooling rates calculated when 
dust and C0 2 are treated simultaneously (labeled “simultaneous 
gas & dust”) 

r. = 0.5 to 1 .0 as they do when dust opacity is changed 
from t v = 0.2 to 0.5 or from 0.0 to 0.2. This occurs because 
line wings are effective in optically, thin media. Once large 
dust opacity is reached, further increases in dust opacity 
do not affect line wing optical depths as much. 

Note that the error increases with dust optical depth, 
even when dust cooling exceeds C0 2 cooling. The expla- 
nation lies in the two-fold nature of dust-C0 2 interaction, 
and underlies the importance of its inclusion even for 
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Fig. 3. Ratio of the 15pm region (12-18 pm) thermal cooling 
rates computed ignoring dust-C0 2 interaction (labeled “gas 
only + dust only” in Fig. 2) to those computed including dust- 
COj interaction (labeled “simultaneous gas & dust" in Fig. 2) i 
for r,(dust) of 0.2, 0.5 and 1.0 at winter polar latitudes. Hence, ' 
these curves represent the overestimation error of computing 
cooling rates in the 12-18 pm wavelength interval without pro- 
per simultaneous treatment of dust and C0 2 , as a function of I 
altitude and dust opacity 


calculations for dust storm conditions. For low dust opac- 
ities, dust-COj interaction involves the effect of dust in 
choking off emission in the line wings of the C0 2 15 pm 
band. However, for large dust opacities, it is C0 2 that 
is the parasite. For large dust opacities, dust cooling is 
dominant and as much as half of all dust cooling occurs 
in the 1 5 pm wavelength interval (Fig. 1). However, dust is 
unabte to cool effectively in the 1 5 pm wavelength interval 
because emitted photons are quickly reabsorbed by C0 2 . 
In fact, the parasitic effect of C0 2 grows as dust opacity 
increases, as shown in Fig. 3. The error reaches an asymp- 
tote at large dust opacities. 


Summary and conclusions 

As much as half of all atmospheric cooling owing to ther- 
mal emission by dust occurs in the 12-18 pm wavelength 
range. Since this is also the wavelength region where most 
cooling owing to thermal emission by C0 2 occurs, it is 
important to consider the spectral overlap between dust 
and C0 2 . C0 2 infrared cooling and dust infrared cooling 
are both overestimated if C0 2 and dust cooling are com- 
puted separately and then summed. Indeed, ifdust-C0 2 
interaction is ignored, then the total 12-18 pm cooling is 
overestimated by almost 30% for minimal dust loading 
at winter polar latitudes. Furthermore, for one vertical 
optical depth of dust at winter polar latitudes, a simple 
sum of dust cooling and C0 2 cooling from 12-18 pm will 
overestimate the actual total 12-18 pm cooling by over 
40% near the surface. Hence, dust-C0 2 interaction is 
important not only for calculating C0 2 cooling rates, since 
dust opacities overlap the C0 2 1 5 pm band, but also for 
calculating dust cooling rates, since the 15 pm band of 
C0 2 overlaps a significant portion of the infrared region 
where dust cooling is most effective. In other words, even 
when the Martian atmosphere is very dusty, when CO : 
cooling and heating rates are minor compared with dust 
(e.g. Lindner, 1993), the interaction of C0 2 will sig- 
nificantly affect the ability of dust to cool the atmosphere. 
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